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We re-examine the wave function of two-dimensional massless Dirac electron in Q-(BEDT-TTF)2l3 
consisting of four molecules A, A', B and C in a unit cell, using a tight-binding model. We find 
zerolines in the Brillouin zone, on which the component of the wave function becomes zero for B or 
C sites. The zerolines, which are bounded by two Dirac points at ±ko and pass through the M- or 
Y-points, result in a fact that the density of states of the B site exhibits no the Van Hove singularity 
near the energy of the Dirac points. By taking account of the on-site Coulomb interaction within 
the random phase approximation, wc examine the spin fluctuation in order to investigate properties 
of the nuclear magnetic resonance for temperatures T > 50K. In the region for 100 < T < 300K, it 
is shown that the Knight sift for B-site monotonously decreases with decreasing temperature, owing 
to lack of the Van Hove singularity, while it shows a maximum for the other sites (A, A' and C 
sites). In the region for 50 < T < 100K, it is shown that the Knight sift is convex downward and the 
Korringa ratio increases with decreasing temperature for B-site. Such a behavior originates from 
the ferrimagnetic spin fluctuation related to the zerolines. These results are consistent with those 
of the nuclear magnetic resonance experiments. 



I. INTRODUCTION 

Molecular conductor a-(BEDT-TTF) 2 I;r- has highly 
two-dimensional electronic system in the plane of BEDT- 
TTF 1 / 2 molecule owing to layered structure with the 
plane of Ij anion, and has brought much interest by 
the variety of electronic states, such as a charge ordered 
stated, a superconducting state in the presence of charge 
ordering, and a zero gap state (ZGS) with a massless 
Dirac electron 3 -. 

The charge ordered state was suggested theoretically 
using an extended Hubbard models^ 6 -, and was confirmed 
by NMR experiment 7 -. The superconducting state with 
the charge ordering under uniaxial pressure along the 
stacking axis (a-axis)^was also investigated theoretically 
using the extended Hubbard model 9 -. A narrow gap state 
(NGS) was suggested to explain both anomalous increase 
of Hall coefficient and almost-constant resistivity with 
decreasing temperature at high pressures^. From the 
calculation of the NGS, it was obtained that the charge 
gap disappears at high pressure, leading to the density 
of state (DOS) vanishing linearly at the Fermi energy^. 
Furher, the ZGS with a Dirac cone in energy dispersion 
was found theoretically^ using the transfer energies of 
ref. Il3l . and was also confirmed by the first principle 
calculation o 14 ' 15 . The energy spectrum near the Fermi 
energy exhibits two tilted Dirac cones, which are de- 
scribed by the tilted Weyl equatio n 16 ' 17 . The tilt of the 
Dirac cone has been confirmed by a comparison between 
the theoretical and experimental results for the temper- 
ature (T) dependence of the Hall coefficien t 18 ' 19 and the 
angular dependence of the magnetoresistanc o 19 ' 20 . From 
the calculation of the variation of the Dirac point, the 
emergence of a pair of massive Dirac points is predicted 
in the charge ordering state at low pressures^, the merg- 
ing of two massless Dirac electrons is shown at extremely 



high pressures^ 6 -, suggesting robustness of Dirac electrons 
against uniaxial pressures. Effects of electron correla- 
tion for long range Coulomb repulsion in the tilted Dirac 
electron system have been investigated theoretically in 
the absence of magnetic fiel d 22 ' 23 and in the presence of 
magnetic field 2 ^ . As for the short range part of Coulomb 
repulsion, on the other hand, effect of fluctuations in ZGS 
has not yet investigated theoretically in the context of 
inequivalence of BEDT-TTF sites (A = A' ^ B ^ C as 
shown in Fig. 1). 

The property of the wave function plays important 
roles in the DOS and the electron-hole excitation in 
a- (BEDT-TTF) 2I3 . The anomalous behavior of the 
Bloch wave functions exists in the vicinity of the Dirac 
point. The momentum dependences of the velocity ma- 
trix and the charge density exhibit a singularity at the 
point i^ 6 - The angular dependence of the wave functions 
for each site reveals a fact that the absolute value be- 
comes zero for B site or C site at a special direction of the 
point^. The theoretical calculation of Knight shift K a 
and (l/TiT)cr^. has been performed based on a tight- 
binding model for a- (BEDT-TTF) 2 13 where K a oc T and 
(l/TiT) a oc T 2 owing to the Dirac cone spectrum. The 
Knight shift for the respective site reveals a relation, 
Kq > > Kb, which is consistent with experimen- 
tal results^ originates from both the tilt of Dirac cone 
and property of wave function in the vicinity of the Dirac 
point. Thus, the inequivalence of BEDT-TTF sites ob- 
served in NMR experiments reveals inner degree of free- 
dom of Dirac electron in molecular conductors £5r— 

However, the tight binding model is not enough to 
explain the details, e.g, the non-linear temperature de- 
pendences of K a as shown in the following experiment. 
In high temperature region (T > 100K), the Knight 
shift K a with decreasing temperature monotonously de- 
creases for a = B site, while it exhibits a maximum 



for a — A and C sites. In medium temperature region 
(50 < T < 100K), the Knight shift is convex downward 
with decreasing temperature for B site, while the compo- 
nents for A and C sites exhibit linear T-dependencesi 25 i 27 
For a local NMR relaxation rate 1/Ti, the difference in 
(\/TiT) a in a = A, B, and C sites, is small in low tem- 
perature region (T < 50K), i.e., all components are con- 
vex downward with T 2 -dcpcndences for T < 100A'. As 
a result, the Korringa ratio (l/TiTK 2 ) a for a = B site 
clearly increases with decreasing T in medium temper- 
ature region, and exhibits an inequality (1/T\TA 2 )b > 
(1/T 1 TA 2 ) A > (1/TxTA 2 )^ . 

In the present paper, we re-examine the wave func- 
tion in a- (BEDT-TTF) 2 I 3 using a tight-binding model 
where the transfer integrals are given by the first prin- 
ciple calculation-^, and calculate the spin fluctuation 
within the random phase approximation on the on-site 
Coulomb interaction to investigate K a and {\/T\T) a for 
T > 50K. Based on the formulation in §11, we demon- 
strate following new results in §111. We find zerolines in 
the two-dimensional Brillouin zone, where the wave func- 
tion is zero for the components of B or C site. They are 
bounded by the two Dirac points at ±k . The zerolines 
passing through the M- or Y-points give rise to the ab- 
sence of two Van Hove singularities for only the B-site. 
In the high temperature region for 100 < T < 300K, 
the present numerical results on T-dependences of K a , 
are consistent with the experimental results^—, where 
existence (absence) of Van Hove singularities for a = A 
and C sites (a = B site) is essential. In the region for 
50 < T < 100K, it is shown that Ab is convex downward, 
(1/TiTA 2 )b increases with decreasing T, and then the 
inequality (l/T 1 T J Pf 2 ) B > (l/TiTtf 2 ) A > (1/T 1 TA 2 ) C 
is reproduced, owing to the ferrimagnetic spin fluctua- 
tion related to the zerolines and enhanced by the on-site 
Coulomb interaction. These results are also consistent 
with those experimental results 2 ^— . In §IV, summary 
and discussion are given. 



II. FORMULATION 

The model describing the two-dimensional electronic 
system in a-(BEDT-TTF) 2 I 3 is shown in Fig. \£±21%L. 
The unit cell consists of four BEDT-TTF molecules on 
sites A, A', B and C, where A is equivalent to A' so 
that inversion symmetry is preserved, while the sites A, 
B and C are inequivalent. There are six electrons for the 
four molecules in a unit cell, thus the bands are 3/4-filled. 
On the basis of the HOMO orbitals of these sites^^, these 
electrons are described by a Hubbard model with the on- 
site Coulomb interaction U, 

ict 
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FIG. 1: The model describing the electronic system of a- 
(BEDT-TTF) 2 I;r*»22 1 2&. The unit cell consists of four BEDT- 
TTF molecules A, A', B and C with ten transfer energies. The 
a- and &-axis in the conventional notation correspond to the 
y- and x-axis in the present paper. The molecules A and A' 
are equivalent in the presence of the inversion symmetry in 
the ZGS. (Color Online) 

where i,j denote indices of a given unit cell, and a,/3(= 
A, A', B and C) are indices of BEDT-TTF sites in the 
unit cell. In the first term, a\ aa denotes a creation oper- 
ator with spin cr(=t, 4-) and Ua-jp is the transfer energy 
between the (i, a) site and the (J, j3) site. Throughout 
the paper, h and the lattice constant a are taken as unity. 
Hereafter, the energies are given in eV, and the tempera- 
ture is also given by k^T in eV, i. e. leV =10 4 K, where 
is the Boltzmann factor. 
The transfer energies at finite temperature T is esti- 
mated by the interpolation formula 2 ^ 

tx(T) = t x (LT) + (t x (RT)-t x (LT))(T-LT)/(RT-LT). 

(2.2) 

The transfer energies tjf(RT) and tjf(LT) 
are given by the 1st principle calculation^ 
where t al (RT) = -0.0101, i o2 (RT) = 
-0.0476, * a3 (RT) = 0.0093, t bl (RT) = 0.1081, i 62 (RT) = 
0.1109, t 63 (RT) = 0.0551, t 64 (RT) = 0.0151, t aV (RT) = 
0.0088, t a3 '(RT) = 0.0019, t a4 '(RT) = 0.0009, and 
tai(LT) = -0.0267, i a2 (LT) - -0.0511, t a3 (LT) = 
0.0323, t 6 i(LT) = 0.1241, t b2 (LT) = 0.1296, t b3 (LT) = 
0.0513, t 64 (LT) = 0.0512, t al ,(LT) = 0.0119, t a3 '(LT) = 
0.0046, t a4 '(LT) = 0.0060, , where low temperature 
LT = 0.0008 and we put room temperature RT = 0.03. 

The Hamiltonian is diagonalized numerically for a 
given k in each spin subspace, according to 



7jegffo-(k) <^3 7 q(k) = £ 7(T (k) d a -y a (k), (2.3) 

0=1 

where £ 7<T arc the eigenenergies ordered such that, 
Cia(k) > 6 CT (k) > 6a(k) > &a(k) ( 7 = 1,2,3,4 is the 
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band index) , and d aia (k) are the corresponding eigenvec- 
tors. The T-depcndence of the chemical potential owing 
to the electron-hole asymmetry^ is taken into account 
in the present calculation, although it has not been con- 
sidered in [28l 

The bare susceptibility on the site representation is 
given byii 

= -— ^Ga^(k + q,e„ +cj ; )G, 9Q (k,e Il ) (2.4) 

kn 

_T_ V *W»(k,q) ^ , 

n l ^ (i( £n + cj t ) e 7 (k + q))(ie n - e 7 (k|r ' 

where the bare Green function on the site representation 
is defined by 

G aP (k,e n ) = ^d Q7 (k)^ 7 (k)/(ie„ -£ 7 (k)) (2.6) 

7 

with e„ = (2n + l)-7rT and loi = 2nnT are the Matsubara 
frequencies, and the form factor is defined by 

F a p(k, q) = F$ (k, q) + F§> (k, q), (2.7) 
^ } (k,q) = £ ( 2 - 8 ) 

7=1,2 

^?(k,q) = £ O k ><l)' (2-9) 

7=1,2 

i^/(k,q) = d al (k + q)d Pl (k + q)*d 0Y (k)d aY (k)*, 

(2.10) 

with 7^7 indicating another band within the conduc- 
tion (7 = 1) and valence (7 = 2) bands. The form fac- 
tor reflects the character of each site for the intraband 
(F^(k, q)) and interband (F^(k, q)) fluctuations. 

In the presence of spin symmetry, the longitudinal spin 
susceptibility, \ s , and the transverse spin susceptibility, 
X , are given byii 

X s =X ± = (l-X°U)- 1 X° (2.11) 

with U = UI and / is the unit matrix. 

The knight shift K a (the local spin susceptibility) and 
the local NMR relaxation rate (l/Ti) Q with A, A', B and 
C sites are given by^ 

K a - Ex«/^ >°)> (2.12) 

P 

(l/Ti) Q - Tlim Vlm[xl(q,^)]/^. (2.13) 

u— >0 * — ' 

q 

III. RESULTS 

A. Zeroline and density of states 

Figure 2(a), (b), (c) and (d) show absolute values of 
the wave functions of the conduction band (7 = 1) for 



(b) 




FIG. 2: Absolute values of the wave functions of the con- 
duction band (£ = 1) for the components of B site ( 
|dsi(k)dg 1 (k)|(a) and and their coutour plot(b)) and those 
of C site ((c) and (d)) in the Brillouin zone. The zerolines of 
the conduction band are bounded by two Dirac points at ±ko 
as shown in the dark region of (c) and (d), and in schematic 
figure (e) for B site (green solid line) and C site (red dashed 
line). (Color Online) 



the components of B site (Figs 2(a) and (b)) and C site 
(Figs 2(c) and (d)). We find "zerolines" which are curved 
lines where absolute value of the wave function is zero, 
only for the components of B and C sites. The zerolines 
for both B and C sites are bounded by two Dirac points 
at ±ko in the Brillouin zone, where the wave functions 
are discontinuous. In the conduction band, the zeroline 
for B (C) site passes through the M- (F-)point as shown 
in Fig 2(e). In the valence band, on the contrary, the 
zeroline for B (C) site passes through the Y- (M-)point. 
The zeroline for B (C) site in conduction band coincides 
with the zeroline for C (B) site in valence band in the 
vicinity of the Dirac points. , while these are slightly 
different far from the Dirac points, since the conduction 
and valence bands are not equivalent in a-type organic 
conductors. 

In Fig 3(a), the conduction band and the valence band 
are shown The chemical potential with u = is situated 
on the Dirac points at T = 0, while it decreases with 
increasing temperature owing to the asymmetry of the 
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FIG. 3: (a) Energy dispersions in the Brillouin zone for the 
conduction and valence bands (7 = 1 and 2), which has a 
contact point with the tilted Dirac cones at ±ko. The chem- 
ical potential is given by with to — 0. There are the saddle 
points at the M-point in the conduction band and at Y-point 
in the valence band near u = 0. (b) The density of states 
p a near the energy of the Dirac point (cj — 0) at T — for 
the components of a = A site (blue dotted line), a — B site 
(green solid line), and a = C site (red dashed line). (Color 
Online) 



B. Knight shift and NMR relaxation rate 

Figure 4(a) shows temperature dependences of the 
Knight shift K a for a = A, B, and C sites with 7 = 0.12 
and K° with U = 0. The value of 7 = 0.12 is cho- 
sen so that the temperature dependences of K a in the 
present calculation reproduce those of experimental re- 
sults for T > 50K. The Knight shift K° (with (7 = 0) 
for T < 0.015 has been calculated in ref. HH, which 
gives slightly different results due to ignoring the In- 
dependence of the chemical potential. In high tempera- 
ture region (100 < T < 300K), With decreasing tempera- 
ture, the Knight shift monotonously decreases for S-sitc 
while there is a maximum in the temperature dependence 
of the Knight shift for A and C sites. Such behavior is 
independent of U originates from the combined effect of 
the Van Hove singularity and the zerolines as discussed in 
Fig 3. In medium temperature region (50 < T < 100K), 
on the other hand, it is shown that the Knight shift is 
convex downward with decreasing temperature for B-site 
with U = 0.12, while the components for A and C sites 
exhibit linear T-dependences in this region. Such a effect 
of electron correlation originates from the ferrimagnetic 
spin fluctuation described later. Figure 4(b) shows tem- 
perature dependences of (l/T\T) a for a = A, B, and 
C sites with 7 = 0.12 and (l/TiT)° with 7 = 0. It is 
shown that (l/T\T) a exhibit linear T 2 -dcpcndcnccs with 
T < 100/v approximately, and all components are weakly 
enhanced by U. Figure 4(c) shows temperature depen- 
dences of the Korringa ratio (l/TiTK 2 ) a for a = A, B, 
and C sites with U = 0.12 and (l/TiTX 2 )° with (7 = 0. 
It is clearly shown that the Korringa ratio for B site is 
strongly enhanced by U and increases with decreasing T, 
owing to the anomalous T-dcpcndence of Kb- The elec- 
tron correlation effect results in the enhancement of the 
following inequality, 

(1/7iTX 2 )b > (l/TiTX 2 ) A > (l/TiTK 2 ) c . (3.1) 



C. Ferrimagnetic fluctuation in spin susceptibility 



conduction and valence bands*^ The saddle points close 
to the Fermi energy are seen at the M-point and at Y- 
point, respectively. The saddle points give rise to the 
Van Hove singularities w — 0.012cV in the conduction 
band, and at lo = — 0.026eV in the valence band. Such a 
singularity appears in the components of density of state 
at A and C sites as shown in Fig 3(b), where the energy 
of the Dirac points is u> = 0. For the component of B 
site, however, these two Van Hove singularities disappear, 
since the zeroline for B site passes through the location 
of the saddle point. 



We examine static spin susceptibility at low energy by 
choosing T = 0.02 and oj = 0. In Fig 5(a), the mo- 
mentum dependences of the diagonal components of the 
spin susceptibilities Xa a an d the bare susceptibilities x° a a 
(U = 0) for a = A, B, and C are shown for a = A, B, 
and C with U = 0.12. The diagonal components are 
positive in the Brillouin zone and are enhanced by U. 
The momentum dependences of the off-diagonal compo- 
nents of the spin susceptibilities Xap f° r ( a > P) = (A B), 
(B, C), and (A, C) with U = 0.12, and the bare suscepti- 
bilities X % for (a,/3) = (A, B), (B,C), and (A, C) with 
7 = are shown in Fig 5(b). It is found that Xbc an< ^ 
Xab are negative in the Brillouin zone, while Xac has 
both positive and negative value. Those absolute val- 
ues are enhanced by 7. Especially, Xbc P^ a y s the most 
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FIG. 4: (a) Temperature dependences of the Knight shift K a 
with U = 0.12 (thick line) and K° (thin line) for a = A (red 
dashed line), B (green solid line), and C sites (blue dotted 
line), respectively, (b) Temperature dependences of (l/TiT) Q 
with U = 0.12 (thick line) and (l/TiT)° (thin line) for a = 
A (red dashed line), B (green solid line), and C sites (blue 
dotted line), respectively, (c) Temperature dependences of 
Korringa ratio (l/TiTTf 2 ^ with U = 0.12 (thick line) and 
(l/TiTA' 2 )° with [7 = (thin line) for a = A (red dashed 
line), B (green solid line), and C sites (blue dotted line), 
respectively. (Color Online) 



important role for the anomalous T-dcpendences of Kb, 
since IxbcI at q = is strongly enhanced by U, owing to 
product of negative Xbc w ^ n positive diagonal terms in 
the RPA process such as 

[x s hc = [(l-X°Ur 1 X°}BC 
= Xbc + ^XbbXbc + ^XbcXcc + ' ■ ' ■ ( 3 - 2 ) 

Negative off-diagonal susceptibilities, Xbc ana - Xab' 
and positive off-diagonal susceptibility, xlc ; at q = in- 
dicate the fcrrimagnctic spin fluctuation where the spin 
on B site tends to be opposite to that of the other sites 
as shown in Fig 5(c). The Knight shift is a sum of the 
diagonal and off-diagonal spin susceptibilities, where the 
former are dominant at high temperatures since the spin 
fluctuation exhibits local character. The latter become 
important with decreasing temperature, resulting in the 
anomalous T-dcpendcnce of Kb which is convex down- 
ward at low temperatures. 

The negative value of Xbc at q = originats from the 
form factor for the intcrband fluctuation, F BC ! (k, 0), since 

the form factor for intraband fluctuation, F^,(k, 0), is 
positive for any k. The momentum k dependence of 

(2) 

F^(k, 0) in the Brillouin zone is shown in Figs 6(a) and 

(2) 

6(b). Although F^(k,0) can take positive or negative 

(2) 

values, the numerical result show F^{k, 0) is negative 
or zero for any k. It is not a self-evident result owing 
to the phase structure of the wave function as discussed 
later. The white curved lines in Fig. 6(b) correspond to 

(2) 

the zerolinc for B and C sites, where F BC !(k, 0) = 0. 

(2) 

In order to show the reason of .F B( j(k, 0) < 0, we 
investigate the phase structure of the wave function. 
Properties of wave functions in the vicinity of the Dirac 
points have been investigated by analying in terms of the 
Luttingcr-Kohn representation 2 ^. The analysis of wave 
functions in the present paper, on the other hand, is 
based on the Bloch representation. Figure 7 shows nu- 
merical results of wave functions for B and C sites and 
for conduction and valence bands, dBi(ko')*i ^B2(ko')j 
^Ci(ko') and dc2(ko')"\ on the Gauss plane, here |ko' — 
k | = 0.057T, i.e., the momentum k ' circles the Dirac 
point k . In order to avoid arbitrary phase factor of the 
wave functions for each band 7 and each momentum ko , 
we choose dA-yO^o') is real and positive. Such a represen- 
tation for wave functions was used in ref. 3% to examine 
the rotation of the base, d Qj i and d a ,2, around ko. The 
inset shows ^-dependences of absolute values of the wave 
functions for the conduction band. The similar result 
was shown in|28| although there is slight difference in the 
choice of the horizontal axis. Based on this numerical 
results, the wave functions in the vicinity of the Dirac 
points, arc approximately given by 

d m (6) = e w '\d B i\[e i{e+e ^ ) - e W™+ e ^], 
d cl (6) = e w i\d cl \[e^+ 9 ^ -e^+Oci)], 
d B 2(0) = e ie2 \d B 2\[e i{e+eB2) -e^+^l 
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d C2 {9) = e Wl \d C 2\[e i{e+9c2) - e ^ BZ+ec2 \ (3 3) 

where 9 is the angle of the vector ko' — ko measured from 
/c^-axis. The values of wave functions at 6 = are deter- 
mined by #bi- $ci; $B2 and 9q 2 ■ The angle of the zeroline 
for B (C) site in the conduction band is 9bz (#cz), where 
kfei(#)| = (\dci{0)\ =0) as shown in the inset of Fig. 
7, while the zeroline for B (C) site in the valence band 
corresponds to 9cz (&bz)- The arbitrary phases of the 
wave functions for each band and each momentum, 9\ (9) 
and #2(0), disappear in the form factor. Using above 
wave functions with the parameters given by the numer- 
ical calculation, we obtain F^(k°,0) < at arbitrary 
9. 

If #cz = #bz +7r (close to the present numerical result) , 

Fj^(k°,0) = -8cos<9 sin 2 ((9 - BZ ) (3.4) 

withflo = -0bi+0 C i-0 C 2+0b2- We obtain F^(k°,0) < 
in the wide region of parameter, — 7r/2 < #0 < tt/2. 

Thus the condition of Fg^(k°, 0) < is robust when the 
zcrolines of B and C sites extend in opposite directions 
each other. When 6>cz = #bz, on the other hand, 

F^(k°,0) = 8cos6> (1-cos((9-<9bz)) 2 - (3.5) 

We obtain F B 2 (~, (k° , 0) < with 9q m the present numer- 
ical result. Thus the angle between the zerolines of B 
and C sites is an important factor for the ferrimagnetic 
fluctuation in the Dirac electrons of a-(BEDT-TTF) 2 I 3 . 

IV. SUMMARY AND DISCUSSION 

In summary, we examined the wave function and the 
spin fluctuation in o>(BEDT-TTF) 2 l3, using a tight- 
binding model and the on-site Coulomb interaction 
treated within the random phase approximation. The 
effect of electron correlation on K a and {\/T\T) a with 
A, A', B and C sites for T > 50K was investigated paying 
attention to the inequivalence of these sites (A = A' ^ 
B^C). 

We found that zerolines, where the wave function is 
zero for the components of B or C site. They give the 
vanishing of two Van Hove singularities near the energy 
of the Dirac points only for S-site component of DOS. 
Existence (absence) of the Van Hove singularities plays 
essential role for the T-dependences of K a . For the 
high temperature region of 100 < T < 300K, with de- 
creasing T, Kb decreases monotonously while Aa and 
Kc exhibit a maximum. In the region for 50 < T < 
100K, Kb is convex downward, and (1/TiTA" 2 )b in- 
creases leading to an inequality of the Korringa ratio, 
(I/TiTA^b > {\/T\TK 2 )x > (l/TxTA^c, was ob- 
taincd. These results are consistent with those of exper- 
iment for T > 50K^M . 



It is found that the anomalous T-dependence of Ab 
is ascribed to the ferrimagnetic spin fluctuation which 
is enhanced by the on-site Coulomb interaction. Such 
a fluctuation describes a spin on B site being opposite 
to the other spins on A, A' and C sites. The ferri- 
magnetic spin fluctuation originates from the interband 
fluctuation mainly between B and C sites. The inter- 
band fluctuation relates to the zerolines bounded by two 
Dirac points in the Brillouin zone. Such zeroline does 
not exist in the wave function of graphene, since carbon 
atoms in two sublattice are equivalent owing to the in- 
version symmetry, which corresponds to A and A' sites 
in a-(BEDT-TTF) 2 I 3 . Thus the present results reveal 
that the inequivalence of BEDT-TTF sites play impor- 
tant roles for observables in NMR as an inner degree of 
freedom of Dirac electron in molecular conductor with 
the short range Coulomb interaction. 

There remain problems to be clarified. In the low 
temperature region for T < 50, K a for all compo- 
nents are convex downward with decreasing T—— , and 
(l/TxT) Q exhibits complex T-dependence for very low 
temperature s 27 ' 31 in the presence of magnetic field per- 
pendicular to the conducting plane. Those behavior can 
not be explained within the random phase approximation 
on the on-site Coulomb interaction. It indicates an im- 
portance of a higher order correction such as self-energy 
correction with the long range Coulomb interaction, since 
at low temperatures, low energy phenomena in the vicin- 
ity of Dirac point is dominant and then scale of length 
is much longer than lattice constant. Finally, we note 
that the valley splitting owing to the pseudo-spin XY 
ferromagnetism in N = Landau states^ may also play 
significant roles at very low temperatures in the presence 
of magnetic field perpendicular to the conducting plane. 
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FIG. 5: (a) Momentum dependences of the diagonal compo- 
nents of the spin susceptibilities Xaa with U = 0.12 (thick 
line) and the bare susceptibilities \aB (thin line) for a = A 
(red dashed line), B (green solid line), and C sites (blue dot- 
ted line), where T = 0.02 and u> = 0. (b) Momentum depen- 
dences of the off-diagonal components of the spin susceptibil- 
ities Xa/3 U = 0.12 (thick line) and the bare susceptibilities 
X% (thin line) for (a, /3) = (A, B) (gold dashed line), (B,C) 
(wine red solid line), and (A, C) sites (purple dotted line), 
where T = 0.02 and ui = 0. (c) Schematic figure for pattern 
of the ferrimagnetic spin fluctuation. (Color Online) 
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FIG. 6: (a) Momentum k dependence of the form factor 
FgQ(k, 0) for the interband fluctuation in the Brillouin zone. 
The contour plot corresponding to (a) is shown in (b) (Color 
Online) 




FIG. 7: Numerical results of Wave functions normalized by 
dA7/|rfA7i| (7 =1)2) on the Gauss plane, for dBi(ko')* (the 
green thick solid circle), dB2(ko') (the green thin solid circle), 
rfci(ko') (the red thick dashed circle) and dc2(ko')* (the red 
thin dashed circle), where |ko' — ko | = 0.05-7T and the angle 
of ko' — ko vector from fc^-axis, 9, moves from zero to 2tt. 
The directions where the values of wave functions move with 
increasing 9 are represented by the arrows on the circles, and 
the values of wave functions at 9 = are indicated by the 
filled or open small circles on the large circles. The inset 
shows ^-dependences of absolute values of the wave functions 
for the conduction band, |dAi(ko')| (the blue dotted line), 
|dBi(ko')| (the green solid line) and |<fci(ko')| (the red dashed 
line. (Color Online) 



